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1. Introduction 

Piecewise linear maps with a single discontinuity have been studied from many points of 
view. Renyi introduced the (3— transformations (maps of the form /3x (mod 1), > 1) 
to study invariant measures and properties of Diophantine approximations, and the idea 
was taken up by Parry amongst others pQ. This led to a body of work on arithmetic 
dynamics where arithmetic properties of the slopes (e.g. being Pisot numbers) lead to 
special cases or transition points [2]. Other approaches to these maps have been recently 
developed by Dajani et al. [3] and Gora [HE]. During the 1980s, piecewise monotonic 
maps with a single discontinuity were used to model global bifurcations, leading to a 
greater understanding of the dynamics through kneading theory and renormalization 
(or induced maps) [6j [7] . In the last few years there has been a renewed interest due 
to applications in non-smooth bifurcation theory. Avrutin et al. have described many 
features of the piecewise linear cases from this point of view [H [9], [10] and maps with 
constant slope have been derived as models of nonsmooth bifurcations of flows [TT1 [T2]. 

In this paper we develop a theoretical approach to understand an observation made 
in |13| . The maps studied are piecewise linear with one increasing branch and one 
decreasing branch. For fixed values of these slopes there is a family of maps having an 
invariant interval and parametrized by the difference in the values of the maps at the 
discontinuity. The numerical experiments of [13] suggest that there are values of the 
slopes such that the Lyapunov exponent (and, as we will show, the topological entropy) 
of the map is constant over a range of choices for the values at the discontinuity. The 
reason for this turns out to be a surprising robustness in the structure of the associated 
invariant measures, but there are still open questions about precisely which slopes admit 
such plateaus. The answer appears to be arithmetic, involving either Pisot numbers or, 
at the very least, more general algebraic integers. 

Before describing the results in detail it is worth noting that it is quite easy 
to construct families of maps with constant Lyapunov exponents and/or constant 
topological entropy, and so we will mention a few of these constructions here to contrast 
with the explanation for the phenomenon in this family. Our aim in the next few 
paragraphs is to show that the results presented here are truly unexpected and cannot 
be understood by trivial mechanisms. 

The most simple way to ensure that entropy and Lyapunov exponents are constant 
in a family is that all elements of the family are smoothly conjugate to some map, 
F : / — > I, say. That is, there is a continuous family of diffeomorphisms : I — > I, 
parametrized by /i, and the family of maps considered is simply G^ = q^ 1 o F o g M . The 
smooth conjugation of all elements in the family G M to F means that provided F has 
positive entropy and well-defined Lyapunov exponent, then all the functions also 
have the same values of these quantities. 

The /3— transformations themselves provide another example. Since the slope 
is always (5 (where defined) the Lyapunov exponent is trivially log/3, and so all 
0— transformations with the same value of /3 have the same Lyapunov exponent. Once 
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again, our examples do not fall into this category, though there is a sense in which our 
maps are closely related to the family of maps 



with s > 1. Provided the parameters are chosen so that an interval is mapped to itself 
then the entropy and the Lyapunov exponent equal logs. Milnor and Thurston prove 
([6], see also Glendinning [13] for an explicit account in this case) that if a piecewise 
monotonic map / with a single discontinuity and one branch increasing with the other 
branch decreasing has positive entropy then it is semi-conjugate to one of the maps P^ s 
and has entropy log s, i.e. there is a monotonic continuous but not necessarily invertible 
map q such that P^ s o q = q o /. Since all the maps we consider have the properties of 
such / they are each semi-conjugate to P^ s for some // and s. The surprise is that over 
ranges of parameters the families are semi-conjugate to maps with the same value of s 
even though the dynamics is changing, i.e. the values of \x is not constant. 

One further way of creating families with constant entropy is a mechanism seen in 
the logistic (quadratic) map. In this case the non-wandering set, at certain values of the 
parameters, can be decomposed into two parts: a fixed part, on which the dynamics has 
positive entropy, and a second component whose dynamics changes with the parameters 
but whose entropy is less than that of the fixed part, which is therefore the entropy 
of the map for all the relevant parameter values. This is the case (for example) in 
the period three window, and more generally when induced maps can be defined for 
higher iterates of the map. In the period three window there are three sub-intervals 
(one containing the turning point) that are permuted by the map and mapped into 
themselves under the third iterate. Between these three intervals there lies an invariant 
set with entropy log 1+ 2 , and the maximum entropy of the third iterate is | log 2, which 
is less than the entropy of the constant part of the non-wandering sets. Glendinning and 
Hall [15] explore this, and the related kneading theory, for piecewise increasing maps. 
A description of topological entropy for these maps is given in |16j . 

The point we are making is that whilst there are a number of simple ways of finding 
families of maps with constant entropy and Lyapunov exponent, the example presented 
here does not fit into these models. Indeed, another way of looking at our results is 
that we have a two-parameter family of maps, and that we are interested in contours of 
constant topological entropy. If these contours were in standard position there would 
be no reason to expect there to be a simple choice of parameters such that a part of a 
contour is a line parallel to a parameter axis. However, in our case this does happen, 
and we believe the reasons reflect arithmetic resonances in the equations. 

In this paper, we restrict ourselves to showing that the phenomenon described 
really does happen. This will involve constructing invariant measures for the maps and 
using these to describe the Lyapunov exponents. As such, we provide further explicit 
examples of the type described more generally by Gora |5]. We will also use kneading 
theory to prove the constancy of the topological entropy of the maps. 




(1) 
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In [T3] we investigated the properties of the map introduced by Varley, Gradwell 
and Hassell (VGH map) in [17] for the study of insect populations. The map has the 
following form 

J ry n , y < c , 
{ ry v n b , y>c. 

The map ([2]) has an exact linearization in terms of the new parameters 

z = 21og(y)/log(r), C = 21og(c)/log(r), (3) 
with y ^ and r > 1, that turns it into a piecewise linear map 

{z n -\- 2, z n <J , , 

(1-6K + 2, z n >£. [V 

A further change of variable (x = z — £) leaves the map in the form 

-sx n + /3, x n > 0, 

with (3 = 2 — 6^ and s = 6 — 1. 

As we showed in [13] the map T(x) is chaotic for s > 1. In the rest of the paper 
we consider only this chaotic regime and restrict our attention to the invariant set 
which is [(1 - s)(3,(3] if (3 > 2, and \J3 - 2s, 2] if (3 < 2. In figure Q] the bifurcation 
diagram as a function of parameter (3 and the corresponding Lyapunov exponent are 
shown for the particular case s — 2. The dynamics is chaotic with positive Lyapunov 
exponent. Interestingly, the numerically computed Lyapunov exponent shows a plateau 
in the parameter range 2 < (3 < 5 on which the Lyapunov exponent appears to be 
constant. The goal of this paper is to show that this is a real phenomenon and not 
a mere numerical artifact. To do so we calculate the Lyapunov exponent analytically 
which in this case reduces to finding a solution of the Perron- Frobenius equation for the 
invariant density. Moreover, we also prove that there is a corresponding plateau for the 
topological entropy; numerical experiments suggest that this coincides with the plateau 
of Lyapunov exponents. 

By kneading theory arguments it is possible to show that the dynamics of the 
map change across the plateau. And also, that none of the simple mechanisms causing 
Lyapunov exponent plateaus that were described above apply in our case. 

The structure of the paper is as follows. In Section [2] we describe the Lyapunov 
exponent plateaus. In Section [3] we explain how to calculate them analytically via 
the Perron-Frobenius equation. We illustrate this method with two simple examples 



and leave the complete proofs to Appendix A and Appendix B Similarly, in Section 



H] we show that the Lyapunov exponent plateaus are accompanied by plateaus of the 
topological entropy and we also prove these making use of kneading theory. In Section 
[5] we summarize our results and propose a list of open challenges to stimulate further 
work. 
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Figure 1. Bifurcation diagram (bottom panel) and numerical estimation of the 
Lyapunov exponent (upper panel) as a function of /? for s — 2. The Lyapunov exponent 
appears to be constant in the range 2^/3^5. 
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2. Lyapunov exponent 

The Lyapunov exponent measures the average exponential rate of divergence of two 
initially close orbits of a dynamical system, with a positive exponent indicating chaotic 
behaviour. In the case of a discrete map /, the Lyapunov exponent of x is 

A(x 0) /)= lim |i^ln|f(x fc )|l (6) 
I k=o ) 

provided the limit exists. 

If fi is an invariant measure for / then the Lyapunov exponent of / with respect to 

H is 

A„ = y ln|/'(x)|d/i(x), (7) 

and for //-almost all points this equals (jSJ). For our map, In = if x < so (J6]) 

becomes 



k=0 



where 9(x) stands for the Heaviside step function. Similarly, equation j7]) becomes 

POO 

= In \s\ / dji{x). (9) 



It is clear from expressions ([5} and (Q that what determines the exact value of the 
Lyapunov exponent is the fraction of time (or number of iterations) an orbit spends, 
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on average, in the region x > 0. Therefore, proving the constancy of the Lyapunov 
exponent in the plateau reduces to proving the constancy of /x([0, oo]). 

Figure [2] shows different plateaus for different values of s. Apparently, only some 
special values of s show a plateau of constant value of the Lyapunov exponent. In 
fact, as figure [3] illustrates, this plateau disappears if we vary the value of s slightly. 
When s6N the Lyapunov exponent remains constant in a range of of varying length 
(depending on the specific value of s) but starting in all cases at (3 = 2. When s takes 
other algebraic integer values, like Pisot numbers (including the golden ratio, the silver 
ratio, etc.), plateaus of constant Lyapunov exponent can also be observed (figure HJ). 
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Figure 2. Lyapunov exponents as a function of f3 for different values of s € N and 
s = (j) = 1±&, 
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Figure 3. Lyapunov exponents as a function of /3 for different values of s close to the 
integer value s = 2. 
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Figure 4. Lyapunov exponents as a function of /3 for the case of s taking the value of 
different quadratic Pisot numbers. 



We have been able to prove that these plateaus are not a numerical artifact by 
determining analytically the Lyapunov exponent for two particular values of s (s = 2 
and s = (f) = 1+ 2 )■ Taking into account fl9]) the determination of the Lyapunov 
exponent is reduced to finding the natural invariant measure p(x) of the attractor. To 
do this, it is necessary to solve the Perron- Frobenius equation for the invariant density. 
This process is explained in the next section. 



3. Lyapunov exponents and the Perron-Frobenius equation 

In general, note that if an invariant measure \x exists with a density p[x) such that 

pb rb 

n([a, b]) = / dp(x) = / p(x)dx, (10) 



then this density satisfies the Perron-Frobenius equation 

_p(y) 

{y\T(y)=x} 

This is a functional equation for the density and there is no standard procedure to solve 
it. However, for some parameter ranges we have been able to pose an ansatz for p(x) 
that we can prove satisfies the Perron-Frobenius equation and is therefore a solution. In 
particular, we find a partition such that by assuming the invariant density is piecewise 
constant, we are able to turn de Perron-Frobenius equation into a simple system of 
algebraic equations. 
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For the case /3 > 2 the Perron-Frobenius equation can be written as 

ip(i(/3 - x)) if x G ((1 - s)/3, (1 - s)(3 + 2), 

p(x)={ p( x -2) + ± P a(J3-x)) if xG ((l-s)/3 + 2,2), (12) 



V(i(/?-x)) if are (2,0). 



We will now illustrate this general method by solving two simple examples for 
particular values of the parameters. 

3.1. Two simple examples for s = 2. 

3.1.1. — 2. If = 2 then the map is continuous and one of the regions in (fl2l) is 
trivial, so 



of ^_J !p(K 2 -*)) ifxG(-2,0), 

/( J S p{x-2) + \p{\{2-x)) ifxG(0,2) Udj 
Considering now the preimages of the intervals I\ = (—2, 0) and I2 = (0, 2) 

T- 1 (/ 1 ) = (l,2), T- 1 (/ 2 ) = (-2,0)U(0,l), (14) 

and assuming the invariant density is piecewise constant as an ansatz, 

f A if x G (-2,0), 
P(X)= 5 if* 6 (0,2), (15) 



the Perron-Frobenius equation ( 1131) gives 

A = is, 5 = ^ + ^5. (16) 
It is easy to verify the simple piecewise constant solution 

P{ ) \ I ifxG(0,2), (17j 

where the normalization has been chosen so that J 2 2 p(x)dx = 2(A + B) = 1. Hence, 
the invariant measure for x > reduces to 

»2 



//([(). 2]) / \dx= 2 -, (U 



and ([9]) gives 



X = hn2. (19) 



(21) 



(22) 
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3.1.2. (3 = 4. The previous case was so easy that the solution could be presented with 
no explanation at all. The case = 4 is a little more complicated and serves to illustrate 
how the complexity of the calculation increases with more complicated dynamics. The 
Perron- Frobenius equation (TP2"]) becomes 

f |p(|(4-x)) if xe(-4,-2), 

p(x)=l p( z -2) + |p(|(4-x)) if x G (-2,2), (20) 
[ \p(\(A-x)) if 16(2,4). 

We consider in this case the following intervals 

h =(-4,-2), 

h =(-2,0), 

h =(0,2), 

/ 4 =(2,4), 

whose preimages are 

T-\h) =(3,4), 
T-\I 2 ) =(-4,-2)U(2,3), 
T-^Ja) = (-2,0) U (1,2), 
T-\h) =(0,1). 

This suggests that we solve ( 1201) by posing the piecewise constant solution 

A if x e (-4,-2), 

B if a: G (-2,0), 

C if x 6 (0,2), 

D ifxe(2,4), 

in which case (120|) implies 

A = \D, B = A+\D, C = B + \C, D = \C. (24) 
The normalization requirement J pdx = 1 in this case reads 

2(A + B + C + D) = 1. (25) 
An elementary calculation shows that this is solved by 

A=^, B = \, C=l D = \, (26) 
and so the invariant measure in x > is 

^•^ll^+fM+H' < 27) 

and the Lyapunov exponent is 

A = |ln2, (28) 
which is equal to the Lyapunov exponent of the previous example. 



p{x) 



(23) 
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3.2. Analytical proofs of two plateaus 

Following the same methodology as in the previous examples it is possible to solve 
the Perron- Frobenius equation as a function of 0. This is possible mainly due to the 
robustness of the invariant density p(x) with varying (3. This allows us to find finite 
partitions of phase space on which the value of p(x) in each interval remains constant 
although the size of the intervals can vary with /3. The robustness of p(x) can be 
observed in figure Q] where the different shades of grey indicate different values of the 
invariant density. Between /3 = 2 and /3 = 3 for instance, we can see that the number of 
intervals remains constant and also does the density in each of the intervals. Moreover, 
the borders of the intervals as a function of /3 are straight lines which further facilitates 
the partitioning of phase space and finding an appropriate ansatz for p(x). 
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P 

Figure 5. Lyapunov exponent as a function of f3 for s = 2. The circles represent 
the numerical estimation of the Lyapunov exponent. The solid line stands for the 
analytical Lyapunov exponent calculated from the Perron-Frobenius equation. 




0.38 




2 2.5 3 3.5 4 4.5 5 



P 

Figure 6. Lyapunov exponent as a function of j3 for s = <p. The circles represent 
the numerical estimation of the Lyapunov exponent. The solid line stands for the 
analytical Lyapunov exponent calculated from the Perron-Frobenius equation. 
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Figure 7. Plots of numerically calculated topological entropy as a function of j3 using 
terms up to and including t 100 in ([55)1 for different values ofseN and s — <j) — 1+ 2 V ^ . 



We have proved the existence of plateaus for the cases s = 2 and s = <fi. However, we 
have not been able to prove them in its entirety because the structure of the invariant 
density gets increasingly complex as /3 grows. The results are shown in figure [5] and 



figure El To be specific, the calculations detailed in Appendix A and Appendix B prove 
the following lemma. 

Lemma 1 Consider the maps defined by with (3 > 2. If s = 2 and 2 < < | then 
the Lyapunov exponent is constant and equal to | In 2. If s = (p = 1+ ^ and 20 < (3 < 4 
then the Lyapunov exponent is a constant and equal to In (f). 



4. Topological entropy and kneading theory 

Figures [7] and [H] shows the numerically calculated topological entropy h of the maps 
as a function of j3 for different values of s. They present plateaus in precisely the 
same positions as the Lyapunov exponent plateaus of figure [2] as shown in figure M 
Thus we conjecture that the Lyapunov plateaus are also parameter intervals with 
constant topological entropy. This assumption makes it possible to investigate other 
possible plateaus in more detail, but before presenting the results we will sketch how 
the topological entropy is calculated. 

Milnor and Thurston [B] developed kneading theory for unimodal maps, and noted 
that by considering the kneading invariant of the map as a power series then if the map 
has positive topological entropy this is equal to minus the logarithm of the smallest 
positive zero of the kneading invariant. Moreover, the kneading sequence (seen as 
power series evaluated at this smallest zero of the kneading invariant) then provides 
a semi-conjugacy to a tent map with slopes having constant absolute value equal to the 
reciprocal of the smallest zero. Glendinning and Hall [15J showed how these results could 
be extended to Lorenz maps (piecewise increasing maps with a single discontinuity) and 
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Figure 8. Plots of numerically calculated topological entropy as a function of j3 using 
terms up to and including t 100 in (pJ5| for different values of s close to the integer value 
s = 2. 



Glendinning [13] describes the theory for maps with a single discontinuity, one branch 
of which is increasing and the other decreasing. Proofs of the results described here can 
be found in j6j IT4] . 

Suppose that / : [—a, b] — > [—a, b] with a, b > 0, is a map with a single discontinuity 
at x — 0, continuous and increasing on (—a, 0) and continuous and decreasing on (0, b). 
Then the standard symbolic description (e.g. [IB]) assigns an address to each point 
x 7^ by A{x) = 1 if x > and A(x) = — 1 if x < 0. The itinerary of a point that is 
not a preimage of zero is then just the sequence 

I(x) = A(x)A(f(x))A(f(x))A(f(x)) .... (29) 

Whilst the itinerary is easy to interpret, Milnor and Thurston [|6J observed that the 
same information is contained in a sequence of coordinates that monitor the slope of 
iterates of the map. It is natural to use a plus sign for an increasing slope and a minus 
for decreasing slope, so let 

9 (x) = -A(x), 6 n (x) = -A(f n (x))9 n ^(x), n > 1 (30) 

for points that are not preimages of zero. Thus the sign of 6 n (x) is the sign of the slope 
of f n+1 at x. Moreover, rather than sequences we can work with formal power series 

oo 

P(x,t) = Y J Vk(x)t k (31) 

k=0 

which is called the kneading sequence of x. For points that are preimages of zero, two 
kneading sequences (the lower and upper kneading sequences) can be defined: 

P(x_,t) =limP(y,t), P(x + ,t) =limP(y,t) (32) 

2/f0 ' 2/4,0 

where limits are taken through points that are not preimages of zero. 
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With the standard lexicographical order, a^t k < bj~t k if 
0, — 1, and cij < bj) the sequences P(x,t) are decreasing functions of x and 
by looking at the difference P(x-,t) — P(x + ,t) this function becomes monotonic and 
continuous as a function of x if t — t*, where t* < 1 is the smallest positive zero of 

P(0_,*)-P(0+,t) = 0. (33) 

Such a value of t* always exists if the topological entropy of / is positive and the entropy 
actually equals h = — log t* (see below) . 

The seminal result of Milnor and Thurston [6], which carries over to discontinuous 
maps of the kind considered here [14] . is that if / has positive topological entropy then 
the function 

a(x) = P^)-P{-a + ^) 

is a semi-conjugacy (monotonic and continuous) from / to a piecewise continuous map 
with slopes having modulus 1/t* and the topological entropy of / is h — — logi*. 

Details can be found in lUE], but from the point of view of this paper the important 
point is that the entropy can be calculated by looking at the zeroes of (|33|) by truncating 
the series at order n. Keeping only polynomial terms up to t n , the maximum error is 
yzj-, and since t E [| , 1) (as the entropy of a two branch map of an interval into itself 
cannot be greater than log 2) this gets small as n gets large. This is how figure [7] and [8] 
were calculated. 

Figures [7] and [9] provide strong evidence that the constant Lyapunov exponent 
plateaus described above, and which we have proved exist if the slopes are (1, — s) for 
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Figure 10. Plots of numerically calculated topological entropy as a function of j3 with 
slopes 1.1 and -2.2 infusing terms up to and including t 100 in ([33)) . 



a variety of choices of s (integers, some Pisot numbers) are concurrent with plateaus 
in the topological entropy of the maps. This conjecture provides a way to attempt 
a further investigation into ambiguous 'almost plateaus' observed for the Lyapunov 
exponents in, for example, the case with slopes (1.1,-2.2). Iteration and root-finding 
algorithms make it relatively easy to compute roots of (1331) to high accuracy. The 
entropy obtained for the case for slopes (1.1, —2.2) is shown in figure [101 In this figure 
the possibility of a plateau is apparent, as it was for the Lyapunov exponent. However, 
a close up of parameter values in the apparent plateau demonstrates that there is 
actually a great deal of structure in the variation of the topological entropy, and so 
we believe that this is not a plateau for the topological entropy of the maps. Assuming 
that the conjectured connection between plateaus in the entropy and plateaus in the 
Lyapunov exponent holds, this implies that the Lyapunov exponent plateaus only exist 
at interesting resonances of the form (1, — s) and are not due to resonances associated 
with the ratio of the slopes more generally. 

We will prove the existence of plateaus for the topological entropy in two cases: 
s = 2 and s = <fi. The proof in both cases relies on the algebraic properties of the slope, 
and this should provide a clue to the answer to the more general question about which 
values of s have Lyapunov plateaus. As with the Lyapunov exponents we do not prove 
the existence of the plateau for the largest possible range of values of (3 - to consider the 
entire interval using our techniques would require the separate consideration of more 
and more complicated cases precisely as with the Lyapunov plateaus - but we do prove 
the existence of a plateau. Recall that our numerical results appear to show that the 
Lyapunov plateaus and the plateaus of topological entropy are equal. 

Lemma 2 Consider the map T of |3P with s = 2. If 2 < f3 < 4 then the topological 
entropy of T is log ( 1+ g ) . 
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Proof: We begin by considering the iterates of x = approached from above and 
below respectively, denoting these by + and 0_ respectively If (3 = 2 or (3 = 4 then 
the iterates of 0_ and + are in a finite set which makes it straightforward to use these 
points to create a Markov partition for which the entropy is straightforward to calculate 
(cf. [19J). So assume that 2 < (3 < 4. Then by direct calculation 

T(0+) = /3>0, T(P) = -P<0, T(-f3) = 2 - < 0, T(2 - 0) = 4 - > 0; (35) 
and 

T(0_) = 2>0, T(2)=/3-4<0, T(/3-4) = 0-2 > 0, T(/3-2) = 4-/3 > 0; (36) 
so in particular 

T 4 (0+) = T 4 (0_). (37) 

This makes it possible to compute the difference of the kneading polynomials ( 1531 
explicitly without knowing the full details of each individual series and hence compute 
the zeros. To be specific: using ( 1301) we find 

P + (t) = -l + t + t 2 + t 3 + t 4 P(4:- (3,t) (38) 

and 

P_(t) = l-t-t 2 + t 3 + t 4 P(4- P,t) (39) 

(strictly speaking the unknown polynomial is P((4 — /?)_,£)) and hence using (133]) the 
entropy is minus the logarithm of the smallest positive zero of 1 — t — t 2 . 

□ 

The proof for the plateau with s being the golden mean is similar: 

Lemma 3 Consider the map T of (TJj with s = = (l + V5)/2. 7/20 < (3 < 20+1 then 
the topological entropy ofT is — logt* where t* is the smallest positive zero ofl—t — t 3 . 

Proof: Again, the end-point is easy to deal with separately so assume that 
20 < (3 < 20 + 1 and write /3 = 20 + e, < e < 1. Then by direct calculation 
using the relation 2 = 1 + to simplify nonlinear terms in 0: 

T(0+) = 20 + e > 0, T 2 (0+) = -2 - e(0 - 1) < 0, T 3 (0+) = -e(0 - 1) < 0, 

T 4 (0+) = 2 - e(0 - 1) > 0, T 5 (0+) = 2e > 0; (40) 



and 



so 



Thus 



T(0_) = 2>0, T 2 (0_)=e>0, T 3 ((L) = 20 - e(0 - 1) > 0, 

T 4 (0_) = -2(1 -e) < 0, T 5 (0„) = 2e>0; (41) 

T 5 (0+) = T 5 (0_). (42) 

P + {t) = -1 + t + t 2 + t 3 - t 4 - t 5 P(2e, t) (43) 
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Figure 11. (a) Plot of numerically obtained relation between j3' and (3 for the case 
s = 2. The dashed line corresponds to the linear regression j3' = 0.79116 + 0.60802/3. 
(b) Detail of the fluctuation A = 0' — (c) and (d) Successive zooms of the framed 
areas. 



and 

P (t) = 1 - t + t 2 - t 3 - t 4 - t 5 P(2e, t) (44) 

where strictly speaking we should have checked that 2e is approached from the same 
side in both cases. Hence using ([33]) the entropy is minus the logarithm of the smallest 
positive zero of the difference of the two power series, i.e. of 1 — t — t 3 . 

□ 

In the introduction we claimed that there are some fairly straightforward ways 
of creating families with plateaus of Lyapunov exponents and topological entropy. 
Figure [11] and figure [12] provide evidence for our claim that our results do not fall 
into these 'simple' categories. 

By the kneading theory developed above, any two-branch map such as we are 
considering which has positive entropy is semi-conjugate to a map with slopes having 
fixed modulus equal to 1/t*, where t* is the smallest positive zero of the kneading 
difference ([33]) . Thus each map in the plateau with entropy — logt* is semi-conjugate 
via the function q defined in (134j) to a map of the form 



, 2 + s*x if x < 0, , 

F "M = i ?-,- x ifx>o, (45) 

where s* = 1/t* and we have rescaled (CQ) to match our parametrization for T. Thus for 
each P in the plateau, s* is known and fixed, so there is one parameter to be matched 
and the relationship between (3 and the associated j3' describes reveals the way in which 
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Figure 12. (a) Plot of numerically obtained relation between j3' and (3 for the case 
s = <p. The dashed line corresponds to the linear regression f3' = —0.00512 + 0.88136/3. 
(b) Detail of the fluctuation A = 0' — j3' . (c) and (d) Successive zooms of the framed 
areas. 



the dynamics of the maps in the plateaus change with the parameter f3. Note that 
we believe that in the plateaus the semi-conjugacy is actually a conjugacy, and the 
computer programme which generates figures [11] and [12] provides some evidence for 
this. The diagrams show the relationship between the two parameters j3 for T and j3' 
for F . It is computed by calculating the kneading sequence of + and 0_ for different 
values of /3 in the plateau and then uses a bisection algorithm based on the natural 
order for these sequences [6] [H] for each such /3 to find the value of /3' for which the 
kneading sequence of 0+ equals the known sequence of 0+ for (3. We then check that the 
sequences for 0_ in both cases are the same to good accuracy which provides evidence 
that the two maps have the same kneading series and that there is indeed a topological 
conjugacy (and not just a semi-conjugacy) between the two maps. 

The correspondence between (3 and /?' appears strictly monotonic, indicating that 
no two maps in the plateau have the same dynamics. Moreover, there is a linear scaling 
and the variation about this linear trend certainly appears very complicated if not 
fractal. This suggests that maps in the plateaus have a very rich dynamic structure 
despite having the same Lyapunov exponents and topological entropy. 

5. Discussion 

In this paper we have proved that the plateaus in the Lyapunov exponent observed in 
piecewise linear families of maps in [13] really do exist and have extended the results to 
different slopes (using slopes s = 2 and s — (f) and we have also shown that the plateau is 
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also present in the topological entropy of the maps. We have provided strong numerical 
evidence in figures [TT] and [12] that the dynamics varies in a complicated fashion within 
the plateaus and, in particular, that the dynamics is certainly not constant within the 
plateaus, which would have been an obvious but somewhat trivial way to obtain the 
plateaus. Numerical experiments also suggest that this phenomenon is associated only 
with maps having one branch with slope equal to unity and not with a simple resonance 
(where one slope is an integer multiple of the other for example). Finally we have also 
obtained exact expressions for the Lyapunov exponent at the lower end of the plateaus 
as it becomes different from the constant function on the plateaus. 

The proof of the existence of the plateaus for the Lyapunov exponent uses the 
fact that we are able to calculate invariant densities for the maps explicitly. These are 
piecewise constant on a finite union of intervals similar to those calculated for some 
related examples by Gora [5], and the structure is robust within the plateaus: the 
number of intervals remains constant over several sub-intervals of the plateau, as does 
the value of the density. Only the lengths of the intervals vary. What makes this striking 
is that (as noted earlier) the dynamics changes with changing parameter even though 
the structure of the invariant density does not. The number of intervals on which the 
density changes at discrete values of the parameter, and these appear to accumulate on 
the right hand end point of the plateau. 

The proof of the plateaus for the topological entropy of the maps suggests that one 
of the reasons for the existence of the plateaus is a robust relation between iterates of 
the discontinuity approached from above and from below. This is almost certainly the 
origin of the algebraic properties of those slopes for which we observe plateaus: they 
are Pisot numbers. We conjecture that if a family of maps has a plateau then the slope 
s is a Pisot number. The relations between the iterates of the discontinuity typified 
by f|37|) and fj42|) also change within the plateaus, and we conjecture that these robust 
relationships are related to the robust structure of the invariant densities. 

As well as the conjectures made above there remain a number of questions which 
would help illuminate this phenomenon. 

• Is it true that the plateaus only occur if one of the slopes equals one? If so, can the 
set of slopes which produce plateaus be characterized completely? 

• What characterizes the end-points of the plateaus? Interestingly, it is possible to 
set all the plateaus for the case s £ Nto the same length and in the same parameter 
range by parametrizing the map with a = + |) instead of (3. As figure [T3l 
shows, the Lyapunov plateaus for integer s occur in the range 2 < a < 4. 

• Gora [H |5] provides a general expression for the invariant density of such maps. 
When do these densities have the robust structure described above (what is their 
'bifurcation theory')? Our results also show that the existence of robust structure 
continues to hold below the lower limit of the plateaus, so the robust structure with 
a finite set of intervals is not in itself enough to guarantee the existence of plateaus. 
How are the two cases different? 
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• Do plateaus exist in other piecewise linear models with positive topological entropy 
and a single discontinuity? 

This latter suggestion deserves a little more examination. The piecewise linear models 
described here have one increasing branch and one decreasing branch. There are two 
other possibilities (increasing-increasing and decreasing-decreasing) and in the non- 
increasing case the bifurcation structures are interrelated [20]. However, there is reason 
to believe that neither of the other cases can have plateaus - essentially because initial 
considerations suggest that the only way to have f n (0 + ) = / n (0_) is to have symmetric 
maps, and hence to be in the trivial cases we have effectively ruled out here, where both 
branches have the same absolute value. If this is really the case then the coincidence 
here becomes even more surprising, and lends credence to the suggestion that this 
phenomenon is an example of the interplay of arithmetic with dynamics in similar ways 
to which certain other properties of the general /3— transformations are considered in 
the area of arithmetic dynamics (e.g. [2]). 



■o* s = <f> 




°' 2 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 

a 



Figure 13. Lyapunov exponents as a function of a = ^xi(/3 + §) for different values 
of s £ N and s = 4> = 
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Appendix A. Lyapunov exponent plateau for s = 2 



Appendix A.l. 2 < (3 < § 

The Perron-Frobenius equation ffl2l in this case reads 

( \ p (\(f3-x)) if x e {-0,-0 + 2), 

p(x) = l p ( x -2) + \p(\(f3-x)) if a; G (-0 + 2,2), (A.l) 
[ \p(\((3-x)) if xe (2,/3). 

We start by considering the parameter range 2 < /3 < 5/2. We will do the calculation 
for this case in detail. Consider nine intervals 



h 


= (- 


(3,(3-4), 


h 


= {fi 


-4,-/3 + 2), 


h 


= (- 


(3 + 2,0), 


h 


= (0,0-2), 


h 


= 08 


-2 £) 


h 


= (£ 

\ 2 


,0-1), 


h 


= w 


-1,-/3 + 4), 


h 


= (- 


+ 4,2), 


h 


= (2,0). 


r preimag 


es of these intervals are 


T~ 




= / 9 , 


T~ 


-\h) 


= (0-1,2), 


T~ 


■\h) 


= (-0,-2)U/ 6 , 


T~ 


-\h) 


= (-2,0-4)U(l,§), 


T~ 




= (/3-4,^)U(f,l), 


T- 




= (&±,P- 3) U(j,f), 


T~ 




= (0-3,-0 + 2)U(0 
= (-0 + 2,O)U(£+?,0 


T~ 


-'(Is) 


T~ 




= (0,^)- 



(A.2) 



(A.3) 



2 i) 

z , 2/' 



2), 



We assume the following ansatz for the invariant density 

p(x) = At, if x e hi i = 1, ■ ■ ■ ,9, with A 7 

Considering flA.18j) this gives the equations 

At 
A 2 



A, 



(A.4) 



\A 9 , 



A, 
A, 
A, 
A, 
A 7 
A, 
A, 



\A 7 , 
A x + ±A, 
A x + \A b , 
A 2 + \A b , 

A 2 + i^ 5 , 

A 2 + |A 5 , 

^3 + \M, 



(A.5) 



Lyapunov plateaus 



21 



It is already clear from these equations that A 5 
The additional normalization requirement reads 

{2(3 - 4)A X + (6 - 2/3) A 2 + ((3- 2)A 3 + {/3 

1)A 6 + (5 - 2f3)A 7 + ((3 

+ {13- 2)A 9 = h 

From here it is trivial to obtain the solution 

M =|At, 



A 6 = A 7 = A 8 and A 3 = A A . 



(2-|)A 5 + (f 



2)A 4 
2)A 8 



(A.6) 



A 2 
A 3 
A, 
A 5 
A 9 



34, 
4^4, 

At, 
2 

9' 

A 6 = 
I A, 



(A.7) 



A 7 = As 



2^4, 



which is independent of (3. Finally, the measure of the interval [0, f3] 

M[0, fl) = 0? - 2)A 4 + (2 - ^)A 5 + (| - I) A, 

+ (5 - 2/3)A 7 + (/3 - 2) A s + (f3- 2)A 9 
= ((3- 2)A 4 + (4 - (3) As + {f3- 2)A 9 



[(/?-2) + ~(4-/3) + l(/3-2)L4 4 

2 
3' 



(Ai 



which is also independent of j3 and produces therefore a constant Lyapunov exponent. 

Using the same method it is straightforward to prove that the invariant density is 
the same for the whole range 2 < (3 < 3. Therefore, we can write 



j_ 

18 
1 
6 
2 
9 
1 
3 
1 

V. 9 



if x e (— /5,/3 — 4) 

if x E \fi - 4, -(3 + 2) 

if x e {-(3 + 2^-2) 

if x e {/3-2,2) 

if x e (2,/3) 



(A.9) 



We can check that the measure n{[0,/3]) for this invariant density is again 



^ e( 2,3)([0,/3]) = (/3 



2)^(4-/3)1 + ^-2)1^. 



(A.10) 



By looking at the bifurcation diagram in figure [T] it is possible to see that the invariant 
density changes at (3 = 3. At this particular value the density is 

' i ifxe (-3,-1), 
ifxe (-1,1), 
if xe (1,2), 

if x e (2,3), 



18 
2 
9 
1 

3 
1 
9 



(A.ll) 
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and again ^= 3 ([0,/3]) = § + § + § = §■ 

In the range 3 < < 4 the density is given by 



P/5e(3,4)W 



< - 



J_ 

18 
1 
9 
2 
9 
1 
3 
1 
9 



(A.12) 



if x e (-0,-/9 + 2), 
if i6 (-/3 + 2,/3 - 4), 
ifxe (/5 - 4,/3 - 2), 
if x e (0-2,2), 
if x e (2,(3). 

The point = 4 is again an special point where the density changes. At this 
particular point the density is given by 

if x e (-4,-2), 
if x e (-2,0), 
if x e (0,2), 
if x e (2,4). 



(A.13) 



Finally, the density in the range 4 < < § is 



3e(4,f) 



< £ 



J_ 

27 
J_ 
18 
_5_ 

54 
1 
9 
A_ 

27 
2 
9 

J7_ 

27 

1 

9 

_2_ 

27 



if x e 

if x e 

if x e 

if x e 

if x e 

if x e 

if x e 

if x e 

if x e 



-0,30-16), 
3/3-16,-/3 + 2), 
-0 + 2,3/3-14), 
30-14,-0 + 4), 
-0 + 4,0-4), 
0-4,-0 + 6), 
-0 + 6,2), 
2,-0 + 8), 
-0 + 8,0). 



(A.14) 



It is straightforward to check that the invariant measure /x([0, 0]) for all the densities 
p(x) we have just defined is |. Therefore we have been able to prove the constancy of 



the Lyapunov exponent for the case s 
exponent takes the constant value 



A = - In 2. 

3 



2 in the range 2 < < |. The Lyapunov 



(A.15) 



For 



> | the number of pieces of the invariant density increases considerably 



and its calculation becomes much more involved. In particular, it is very difficult to 
find an appropriate partition on which we can define an ansatz of p(x) that solves the 
Perron- Frobenius equation. 



Appendix A. 2. < 2 

Using the same method it is also possible to study the dependence on of the Lyapunov 
exponent just before the start of the plateau. In this case the Perron-Frobenius equation 
reads 

o(x)-l WW"*)) ifxe (0-4,0-2), 

1 ' ~ s p(x - 2) + Ip(|(0 - x)) if x e (0 - 2, 2), lA ' ibj 
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Making use of the following ansatz 

and taking into account the normalization we get 

if x G (/3 - 4, /3 - 2), 



**Ht» ; f ;-;,_ 2 ; 2 , <^> 

This results in the following Lyapunov exponent for < 2 

A = ^1b2. (A.19) 

Appendix B. Lyapunov exponent plateau for s = = 1+ 2 V ^ 
Appendix B.l. 20 < (3 < 4 

Similarly to the case studied in the previous appendix, the Perron-Frobenius equation 
(fl2"|) is given by 

±p(±(0 - x)) if x E ((3(1 - 0), (3(1 - 0) + 2), 

p(x) = <( p(x - 2) + ±p(±(0 - x)) if x G 09(1 - 0) + 2, 2), (B.l) 
Jp(jCS-x)) if 16(2,13), 

We consider now a partition in eight intervals 

Ji = (0(1-0), 2(0 -20-1)), 

J 2 =(2(0-20-1), 0(1-0) + 2), 

/ 3 = (0(1-0) + 2, 0-20), 

J 4 =(0-20,20-40), 

J 5 =(20 -40, 0(1-0) +4), 

J 6 =(0(1-0) +4, 2), 

J 7 =(2,0(l-0) + 20 2 ), 

J 8 =(0(l-0) + 20 2 ,0). 

The preimages of these intervals are 

T-\h) =(i(00-2),0(l-0) + 20 2 ), 

T-\I Z ) = (0(1 - 0), - 20 - 2) U (2, 1(00 - 2)), 

T-\h) = (0-20 -2, 2(0 -20- 1)) UJ 6 , 

T- 1 ^) = / 2 U (0 + 4(1-0), 0(1-0) +4), 

T- 1 ^) = (0(1 - 0) + 2, 0) U (£=2, + 4(1 - 0)), 

T- 1 ^) =(0-20,^), 

T-^/g) =(0,0-20). 



(B.2) 



(B.3) 
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Assuming the following ansatz for the invariant density 

p(x) = Ai, if x € ij, i — 1, ... ,8, with A 4 = A 5 . 
fIB.ll) gives the equations 



Ax 


= ^8, 




A 2 


= ?^7, 




A, 


= A 1 + 




A, 


= A 1 + 


^6 


A 5 


= A 2 + 




A, 


= ^3 + 




A 7 


= ^5, 




As 


= ^3- 





The additional normalization requirement reads 

(0(1 + 0) - 2(1 + 20))Ax + (4 - /3)(0 + 1)A 2 
+ ((/3 - 2)0 - 2)A 3 + (/3 - 20)A 4 + (4 - /3)(0 + 1)A 
+ (-2 - 0(1 - c/>))A 6 + (/3(1 - 0) + 2<P)A 7 
+ ((/3-2)0-2)A 8 + (/3-2)A> = l. 

From here it is trivial to obtain the solution 

Ai =^A 5 , 
A 2 =jiA 5 , 

4 r _ i 

^ 5 10-40 

A 6 =(j, + i)A 5 , 

A 7 = ^5, 

^8 = ~^A 5 , 

which is independent of /3. Finally, the measure of the interval [0, 0\ is 

which is also independent of /3 and results in the constant Lyapunov exponent 

A = ^ln0 

5-20 Y 
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Appendix B.2. f3 < 20 

To obtain the analytical form of the Lyapunov exponent before the plateau, we consider 
the same Perron- Frobenius equation ( 1B.1I) . However, now we divide the phase space in 
five intervals 

h =(/3(l-0),/3-20), 
h =(/3-20,/3(l-0) + 2), 

h = (/3(l-0) + 2,/3-20 + 2), (B.10) 
h = (£-20 + 2,2), 
h =(2,/3), 

with inverses 

T~\h) =I S , 
T-\h) =h, 

T-\I 3 ) = / 1 U(2(2-0),/3-20 + 2), (B.ll) 
T-\h) = - 20, 0)U(^, 2(2-0)), 
T-\h) =(0,^). 
Using the following ansatz for the invariant density 

p(x)=A i: if x 6 I h i = l,...,5, with A 2 = A 3 (B.12) 

the Perron- Frobenius equation produces the following set of equations 

A 2 = ±A 4 , 

A 3 =A 1 + ±A 3 , (B.13) 
A A = A 2 + ±A 3 , 

Considering the normalization it is straightforward to find 

Ai =j>A 3 , 
A 2 =A 3 , 

As = ^-2(0+0-4) ' ( B ' 14 ) 

Ai =0A 3 , 

The measure of the interval [0, 0\ is therefore 

M[0,/3]) = 2(3-0)A 3 , (B.15) 
and the Lyapunov exponent is 
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